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Abstract 

I Big bang nucleosynthesis (BBN) and the cosmic microwave background (CMB) have a long history to- 
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gether in the standard cosmology. BBN accurately predicts the primordial light element abundances of 



^ . deuterium, helium and lithium. The general concordance between the predicted and observed light element 

abundances provides a direct probe of the universal baryon density. Recent CMB anisotropy measurements, 
' particularly the observations performed by the WMAP satellite, examine this concordance by independently 

<N : 

measuring the cosmic baryon density. Key to this test of concordance is a quantitative understanding of 

CN ■ 

^ . the uncertainties in the BBN light element abundance predictions. These uncertainties are dominated by 



systematic errors in nuclear cross sections, however for helium-4 they are dominated by the uncertainties 



' in the neutron lifetime and Newton's G. We critically analyze the cross section data, producing represen- 

\ tations that describe this data and its uncertainties, taking into account the correlations among data, and 

. . . . 

. explicitly treating the systematic errors between data sets. The procedure transforming these representa- 

' tions into thermal rates and errors is discussed. Using these updated nuclear inputs, we compute the new 

^ ' BBN abundance predictions, and quantitatively examine their concordance with observations. Depending 

\ on what deuterium observations are adopted, one gets the following constraints on the baryon density: 

C§ ! VlYih'^ = 0.0229 ± 0.0013 or Vl^K^ = 0.0216^°;°°^° at 68% confidence, fixing N^^^fj = 3.0. If we instead 

adopt the WMAP baryon density, we find the following deuterium-based constraints on the effective num- 
' ber of neutrinos during BBN: N^eff — 2.78j^Q yg or N^eff = 3.65jl}'3Q at 68% confidence. Concerns over 

^ , systematics in helium and lithium observations limit the confidence constraints based on this data provide. 

BBN theory uncertainties are dominated by the following nuclear reactions: d{d,n)^lle, d{d,p)t, d{p,j)^}ie, 
'^He(a,7)^Be and ^He((i,p)''He. With new nuclear cross section data, light clement abundance observations 
and the ever increasing resolution of the CMB anisotropy, tighter constraints can be placed on nuclear and 
particle astrophysics. 
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I. INTRODUCTION 



The field of cosmology has recently entered a golden age. An age where a global picture of the 
universe is crystallizing because of new precision observations that can test the basic framework of 
the standard cosmological model. With the plethora of new data, it is important to review and test 
the fundamental theoretical pillars of cosmology. These pillars are the theory of general relativity 
and the universal expansion, big bang nucleosynthesis (BBN), and the relic cosmic background 
radiation. 



A. History 

Knowledge of general relativity and the discovery in 1929 by Hubble that the universe was 
possibly expanding , led to the idea that one could extrapolate backwards and conclude that the 
universe was hotter and denser in the past. This idea became what is currently called the "hot 
big bang" model of the universe. Almost 20 years later it was realized that at early enough times, 
the universe would have been hot and dense enough for nuclear fusion to take place. This epoch 
of primordial nucleosynthesis could explain the large abundances of hydrogen and helium seen in 
the universe (and ultimately the trace D, ^He, and ^Li abundances), first explored by Gamow 

a 946), Alpher, Bethe and Gamow (1948), Hayashi (1950), and Alpher, Follin and Herman (1953) 
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The "hot big bang" model also predicted a relic photon background, created when ions recom- 
bined with electrons to form neutral atoms (Alpher & Herman, 1948,1949 0, Q|). In 1965, this 
uniform 3 Kelvin background was detected by Penzias and Wilson for the first time in the mi- 
crowave band 8]. This cosmic microwave background (GMB) offered supporting evidence for the 
"hot big bang" model and stimulated further refinements in the theory of big bang nucleosynthesis 



(Peebles, 1966; Wagoner, Fowler & Hoyle, 1967 |s,|ldl). 

A decade ago, the COBE satellite detected for the first time the 1 : 10^ intrinsic temperature 
fluctuations in the GMB 11]. During the last five years, many more GMB temperature anisotr opy 



measurements have been made (e.g. MAXIMA, BOOMERANG, DASI, GBI, AGBAR 



JJl). The 



atest of these observations being from the WMAP satellite, with its first data release 
in early 2003 Qj. 

These two pillars of cosmology offer a unique probe of early universe physics; while their ul- 
timate concordance depends upon the accuracy of the standard cosmological model and of the 
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observations driving this precision era. These observations are so precise that we can test and 
constrain cosmology in a profound and fundamental way. For reviews of BBN see Schramm and 
Wagoner(1979) 18 1, Yang et al. (1984) Q, Boesgaard and Steigman (1985) Kolb and Turner 



(1990) 



Walker et al. (1991) _22], Sarkar (1996) 



Tytler et al (2000) 



23], Olive, Steigman and Walker (2000) 



25l ] and the Particle Data Group BBN Review by Fields and Sarkar (2002) 



For reviews of CMB theory see White, Scott and Silk (1994) Q, Tegmark (1995) 



30| and Hu and Dodelson (2002) 



231, Van der 
and 



van 




Veen (1998) \2% Kamionkowski and Kosowsky (1999) 
of CMB observations see Wang et al (2002) 32] for a pre-WMAP evaluation and the individual 
group papers mentioned above. 



B. Goals 



Over the past decade, a major thrust of research in BBN has been towards increasing the rigor 
of the analysis. On the theory side, the key innovation was to calculate the errors in the light 
element predictions in a systematic and statistically careful way. This was done u sing Monte Carlo 
analyses (Krauss and Romanelli 1990 [s^; Smith, Kawano and Malaney 1993 [sJli Krauss and 



Nollett and Buries 2000 



Kernan 1995 351; Hata et al. 1996 M: Fiorentini et al. 1998 bl; 

n n 

Cyburt, Fields and Olive 2001 39]; Coc et al. 2002 [40|]), which account for nuclear reaction 



uncertainties and their propagation into uncertainties in the light element abundance predictions. 
These calculations are essential because they allow for a careful statistical comparison of BBN 
theory with observational constraints; in addition, they point the way toward improvements in the 
theory calculation. 

In its standard N^, = 3.0 form, primordial nucleosynthesis is a one parameter theory, depending 
only on the baryon-to-photon ratio rj = n-Q/n^. This is related to the cosmic baryon density; 
assuming that H and "^He are the dominant constituents after the epoch of BBN (pB = ''^h^^h + 
^^He^iHe)! yields the relation: 

2iz.mn,H' . io-„ [1.0 - o.oo7ii86y,l ( ,,,,3 ( J^) ' (D 

where J^b is the current baryon density relative to the critical density, pc = 3H'^ /SttGn. H is 
the current Hubble parameter, usually defined as H = WOh km s^^ Mpc^^ and Gat is Newton's 
gravitational constant. Yp is the primordial, post-BBN mass fraction of baryons in the form of ^He 
and Tyfi is the current temperature of the cosmic microwave background. Since the mass of the 
proton is not the same as the mass per baryon of ^He, Yp appears in eqn. ^ One can see that 
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with the convolution of BBN theory predictions with hght element observations, constraints on the 
baryon density can be placed. The agreement between the various baryon density constraints from 
different light element observations places quantitative limits on their concordance. Deviations 
from concordance, suggests unknown observational, experimental or theoretical systematics. The 
latter possibly indicating the need for new physics in the standard BBN framework. This has been 
extensively explored, the reader is recommended the following incomplete list of reviews 



2^, 
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With little change in observational or experimental data, these bounds have remained relatively 
unchanged over the last few years. 

The recent boon in CMB anisotropy measurements, offers to reshape the cosmological landscape. 
What these observations bring to the table is an independent measure of the cosmic baryon density. 
This independent measurement of the baryon content examines the general concordance of the 
BBN light element abundance theory predictions and their observed values, and tests the basic 
framework of the hot big bang model. It acts as a "tie-breaker" for the various light element 
observation-based baryon density constraints [s? , 
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Key to this test, is an understanding of the dominant uncertainties in the light element predic- 
tions. These uncertainties stem from the systematic errors in nuclear cross sections. We present a 
new procedure for determining cross section representations and their uncertainties and describe 
how they propagate into thermal rates and the light element predictions. With this updated nuclear 
network, we then quantify the concordance between the light element abundance observations and 
their predictions, and the CMB. With this level of concordance set for the standard cosmological 
model, we can test and constrain non-standard models. We use primordial nucleosynthesis and the 
cosmic microwave background together to probe early universe physics spanning times from 1 sec 



to 400,000 yrs after the big bang and beyond. T 
by Cyburt, Fields and Olive (2001,2002,2003) 



lis work follows naturally from the work performed 



3£ 



4e 



47l | and continues with the same 



guard as 



38|. 



the research by Smith, Kawano and Malaney (1993) |3J], and Nollett and Buries (2000) 

This paper is organized as follows: In §^ we describe the formalism of creating representations 
and uncertainties for cross sections and transforming them into thermal rates and uncertainties. 
In § mil we discuss the resulting cross sections and thermal rates and their impact on the light 
element predictions of primordial nucleosynthesis. We then establish the level of concordance 
existing between light element observations, their predictions and the CMB in § IIVI followed by 
conclusions in S IVl 



4 



II. FORMALISM 



In this new age of precision cosmology, it is increasingly important to have an up-to-date and 
accurate theory of primordial nucleosynthesis. Since BBN's uncertainties stem from uncertainties 
in nuclear cross section data, we develop here a rigorous and reproducible procedure for determin- 
ing accurate representations of that data. There are several requirements we wish to impose on 
this analysis. (1.) The representation of the data must be model independent, other than basic 
assumptions of functional form and demanding sufficient smoothness. (2.) The treatment should 
be global, all data is analyzed simultaneously, avoiding operator's discretion and so-called "chi-by- 
eye" systematics. (3.) There should be explicit treatment of a.) the correlations among data in a 
data set and b.) the discrepancies between different data sets' normalizations. These explicit and 
implicit normalization errors dominate over the statistical uncertainties in the data. 

With these goals in mind, we set out to build a framework for representing cross section data 
for the nuclear reactions important for an accurate BBN calculation, seen in tabled To be gm, we 
will discuss the way cross section data is presented, defining notation that will be useful. We then 
present the scheme for determining the best representation of data and the uncertainty in such a 
representation. Finally, presenting the reactions most important for primordial nucleosynthesis, 
and their fits and uncertainties. 

TABLE L Shown in this table are the 12 most important reactions affecting the predictions of the light 
element abundances ("^He, D, '^He, ^Li). 

Reactions 

n-decay 

p{n,j)d 

dip, 7)^He 

d{d,nfKe 

d{d,p)t 

^Ile{n,p)t 

t{d, n)'^He 

3He(d,p)4He 

3He(a,7)^Be 

^Li(p, a^Ue 
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A. Data Sets 



Ideally, a cross section datum contains four numbers: expectation values and uncertainties for 
the cross section and energy. Uncertainties in the energy are typically, negligibly small. One 
difficulty in measuring cross sections is determining their absolute normalization. In addition to 
the statistical error in each point a normalization uncertainty is assigned for a particular data set. 
In many cases this systematic normalization uncertainty dominates over the the statistical. When 
using cross section data, we must take into account the fact that data from a particular data set 
are correlated with each other due to this normalization error. To help with visualization and to 
find the correlation matrix, we define a random variable to draw from to produce a data point: 

Xi = {l + eio)(w + (^iZi)- (2) 

We denote a random variable by underlining it (e.g. x^, z_q, zJ. The random variables, Zjq and z^, 
are assumed to be uncorrelated random variables, with zero mean and unit variance. Notice that 
the mean normalization is unity, we could have allowed another normalization, but opted not to 
because we do not have a reference point to normalize to. In principle, one could use theory to 
determine an experiment's normalization, but we choose the model-independent approach, relying 
on the data as is. Not renormalizing here will lead us to the separate treatment of systematic 
differences between data sets, and the assignment of an overall "theory" normalization uncertainty. 
The expectation values and correlation matrix elements are: 

Exp[xj] = {Xj} = Hi (3) 
Cov\xi,Xj] = (xiXj) - {x^{xj) = (1 + e^)af5ij + e^/ij/Xj (4) 

Generalizing this for multiple data sets we get, 

CinJn = (1 + ^l)(^ijinjn + (5) 

where in denotes the i^^ data point in the n''^ data set. The inverse covariance matrix is: 

2 

(1-1 _ ^injn 

^in,jn n 4-^2^2 '2 „ /„. \2- W 



+ 1+ (t^ 



(1+4) \ CTk^ 

It is this inverse covariance matrix that will be used in the later best fit calculation. In the case 
where the normalization error is small, the covariance matrix reduces to the standard diagonal 
form with the statistical errors as the diagonal elements. In the case where the normalization 
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errors dominate, the inverse matrix becomes: 



.-1 

t'riijn 



(1 + e^>? 



1 



( 





'71 



) 



(7) 



where ag is a typical, albeit small statistical uncertainty. When data sets are large, the second 
term in the parentheses becomes small, thus the covariance matrix again reduces to the standard 
form with slightly inflated statistical errors. Thus, one can see that quite generally, the statistical 
uncertainty is the dominant contribution to the inverse covariance matrix, not the total uncer- 
tainty. When data sets are small, the covariance matrix is highly non-diagonal. However, since 
our prescription combines several data sets, one is often significantly larger than the others, thus 
smaller data sets will have less impact on the fit. 

If we were dealing with one data set, we would not necessarily need this formalism. As noted 
in D'Agostini (1994) 48], it is generally better to treat the normalization error separate from the 
statistical one (e.g. determining a best fit based on the statistical uncertainties alone and adding 
in the normalization error after the fitting process), however we would be ignoring correlations. 
Since we are combining multiple data sets in a meta-analysis, we must include the normalization 
correlations between data points as well as find an effective overall normalization error to add after 
the fitting process. Different data sets may disagree on the shape, or cover different energy regions. 
Therefore, we will continue with the formalism we have laid out. 

B. Creating Representations 

We are interested in determining the best fit parameters for some general, linear parameteriza- 
tion of the data. In each data set, n we have the following data; a position variable (e.g. the energy) 
Xi^^, the expectation value of the function (e.g. the cross section or S'-factor), yi^ measured at Xi„, 
and the covariance between data points Ci^j^. We will assume that each data set is independent 
from all others. 

In the standard treatment, we determine the best fit parameters by minimizing a x^- Foi' 
simplicity, we choose a linear combination of known functions for our parameterization, y{x) = 
X]p ^pXp(x), where Op and Xp(x) are the p^^ of P fitting parameters and fitting functions evaluated 
at X. For example, a polynomial fit (which we will adopt) has Xp(x) = x^. To begin, I will look at 
the case for one data set. We define: 



p 



p 




(8) 



We reiterate here that a caUigraphic C denotes the covariance between data points. 
When determining the best fit parameters by minimizing x^, we can re- write it as: 

p 

p,q=l 

where Op and Cj, ^ are the most hkely values and the inverse covariance between the p*^ and q^^ 
parameters. We note here that an itaUc C denotes the covariance between fitting parameters, not 
the data points. The best fit and its variance are then 

IJi{x) = ap^p{x) (10) 

p 

a\x) = XCp,,Xp(x)X,(a;). (11) 

The most Ukely parameter values are given by: 

p 

q=l 

where 

Cp'i = \Y. ^iJ + ^ii^i)M^3)] ■ (14) 
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Since were are demanding linear fitting functions, the is quadratic in the fitting parameters, 
thus yielding a correlated gaussian probability distribution with the form: 



exp 



4(a — a)^-C ^-{a — a) 

Cid) = — ^ , ^. (15) 

V(27r)^det(C) 



When generalizing this to more than one data set, we have to ask ourselves how do we want 
to weight the data and each data set. If we wanted to rely strictly on the data itself, then the 

is simply the sum of the x^'s from each experiment. This in turn propagates into the fitting 
parameter likelihood distribution as: 

N exp [-1(^-4)^. C(")"'- [a -in) 

'C(a) = TT , F= -■ (16) 

H V(2vr)^det(C(«)) 

This scheme gives more weight to the data sets with more data points, with det(C*^"^) scaling 
like 1/In, where I„ is the number of data points in the n*^ of N data sets. 
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If instead we wanted to treat data sets on an equal footing, then the parameter hkehhood 
distribution takes on the form: 



^ exp 



(a — an 



C{a) = ^y ^ , ^, (17) 



^ti V(2^)^det(C(")) 
Notice the products of hkehhoods has been replaced with a sum of likelihoods in this non-standard 
treatment. We can see with this likelihood, that a analysis becomes more complicated. The 
effective = — 2 In £ is no longer quadratic in the fitting parameters, thus making the distribution 
non-gaussian and possibly multi-peaked. 

Since we are not only determining the magnitude, but the shape of a function, we should rely 
more on the data sets that have more points. Thus the first prescription is appropriate for our 
purposes. The minimum in this prescription is: 

N N 

Xmin = Xmin.n + Yi^L - CLn) ■ C^") (a - a„). (18) 

n=l n=l 

The best fit parameters are still given by equation 1121 but where 

A = E Cr„!,JX,(x,Jy,„+y,„X,(x,J], (19) 

^ N In 

= 2^ Yl C,r;^.JXp(x,JX,(x,J + X,(x,JXp(x,J]. (20) 

n=l i„,j„=l 

Note that if data sets disagree, the minimum per degree of freedom (Xu ~ X^/^) will be 
large, where v is the number of degrees of freedom. With the covariance in the fitting parameters 
depending solely on the covariance among the data, which as discussed earlier depends mainly on 
the statistical uncertainty, the error in the mean, cr(x), is a measure of the statistical uncertainty 
only. When we have a lot of data, this error will be small due to the 1 / \/iV suppression of the error 
in the mean. Thus if we have two data sets with a large quantity of data, but both systematically 
offset from each other, the error will be underestimated. 

This procedure does not take into account the systematic differences between data sets. There 
are various ways of treating uncertainty assignment with discrepant data. The Particle Data Group 
prescription, is to blow up the error in the mean by the factor ^/x^ This has the virtue that it 
does take into account systematic differences and effectively forces the Xu to be unity. Its limitation 
lies in the fact that this scale factor does not cancel out the suppression in the error, thus 

for sufficiently large data sets, this error assignment will still underestimate the true errors when 
using two discrepant data sets. 



When dealing with a one parameter fit or renormahzation where systematics dominate, Cyburt, 
Fields and Olive (2001) showed that the appropriate scale factor for discrepant data is not 
\f~^ js^. This error assignment turns out to be the weighted dispersion about the mean: 

= ^^p^. (21) 

This approach reproduces well the uncertainties when discrepant data are present and has the 
added virtue that it continues to minimize the variance and does not scale with the number of data 
points. Ultimately this error is a measure of agreement between data sets. If agreement is met, 
then the error in the mean will dominate over this error. The limitation of both these methods, is 
that they do not treat the error as if it varies with respect to the position variable x (e.g. energy). 
However, if the differences between data sets is attributable to an unknown normalization error, 
then assuming no energy dependence is appropriate. 
The Nollett and Buries (2000) 3 

compilation does not explicitly calculate systematic uncer- 
tainties. They create samples of mock data, including the intrinsic normalization errors, and adopt 
piecewise, smooth B-spline representations of cross sections, dividing the energy range into smaller 
bins. Each realization is thermally averaged and propagated through the BEN calculation. This 
treatment has the virtue that it has an explicit treatment of the normalization errors, however their 
B-spline fitting procedure does not take into account the correlations between data points. Also, 
this method simply blows up the errors by reducing the number of points contributing to the fit 
in a particular energy bin. This method's main limitation is that it introduces some arbitrariness 
into where the data cuts are placed, and that it is still dealing with a strictly statistical uncertainty 
and not a systematic one. If discrepant data exist such that it lies outside the typical error size, 
then the Nollett & Buries method will tend to underestimate the true uncertainty. In addition, 
the energy correlations of the cross section data are not included (by assumption) in their fitting 
procedure, thus affecting their best fit values. 

It is clear that a procedure is needed to take into account the systematic errors. We will assume 
here, that the systematic errors are purely normalization errors and as such are constant functions 
of energy. We will adopt the Cyburt, Fields and Olive 2001 js^l sample variance as a measure of 
discrepant data. Generalizing its form to take into account the correlations between data, we get: 

In addition, we need to calculate the normalization error inherent to the data. We choose a 
weighting scheme, such that data sets that agree with the fit are given more weight than data 
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sets that disagree, since we have aheady taken into account discrepant data sets. We define the 
intrinsic normahzation error to be: 



^norm — I (^3) 

An 



Here, Xn minimum per datum of data set n, given the best fit parameters. The total 

disc ' norm' 



normahzation error is then the quadrature sum of these two systematic errors 5"^ = Si'j,- + (5^ 



This propagates into our final error as: 

a\x) = (1 + 5'')al,t{x) + S^f^ixf = [(1 + 5^)Cp,, + d^a^a,] Xp{x)X,{x). (24) 

p,q 

There is no unique way to assign a systematic uncertainty. However, any determination being, 
based on the same data, must agree with the overall results of this prescription. This leaves us 
with the question of how can we further improve these uncertainties. There are two ways we can 
improve our errors, (1) we can get new, more accurate and precise data and (2) we can, with 
sufficient reason, exclude data sets, in an effort to remove the cause of the systematic errors. As 
there is not an un-biased way of performing the latter, we rely on the former for the future progress 
of this type of analysis. 

C. Thermal Averaging 

Thermonuclear reaction rates and the reaction networks they belong to, play a key role in nuclear 
astrophysics theory, ranging from stellar interiors, supernovae explosions to big bang nucleosyn- 
thesis. A large base of work has been done in this field. Reaction rate formalism is thoroughly 
reviewed in Clayton's "Principles of Stellar Evolution and Nucleosynthesis" (1983) j^l and Rolf 
and Rodney's "Cauldrons in the Cosmos: Nuclear Astrophysics" (1988) js^]. Compilations of nu- 
clear data and thermonuclear rates began with the pioneering work of William Fowler [sil. Is^. Ib^ . 
A recent update has been provided by the NACRE collaboration |5J]. Recent BBN rate compila- 



3J|, Nollett and Buries 38] and the 



tions have been performed by Smith, Kawano and Malaney 
Cyburt, Fields and Olive ,39] tailored NACRE Q compilation. 

D. Mapping Cross Sections into Thermal Rates 

We want the rate at finite temperature, for 2-body interactions of the type: i + j — > k + I, 
which is Xi+j^k+iiT) = NA{(Ti+j^k+iv) , where A^^i is Avogadro's number. The angle brackets 
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denote thermal averages. We are interested in transforming energy dependent random functions 
into temperature dependent random functions. We define the transformation from one to the other 
as: 

/•oo 

A(T) = / W{E, T)S{E)dE, (25) 
Jo 

where W{E,T) is a weighting function or kernal and is the function we are transforming, 

either the astrophysical S- or R- factor of a cross section. depends on random variables (i.e. 

fitting parameters), and thus is a random function, where the expectation value is Exp[5(-E)] = 
IJ,{E) and the variance is Yar[S_{E)] = a'^(E). 

We want to know how this randomness propagates into A. The expectation value is: 

/•oo 

MA(r) = Exp[A(T)] = / W{E,T)fi{E)dE. (26) 
Jo 

The variance in A is then: 

al{T) = Var[A(r)] = Exp[X{Tf] - Exp[A(T)]2 (27) 

/"OO /"OO 

= / / WiE,T)W{E',T)Cov\S{E),SiE')]dEdE' (28) 
Jo Jo 

/•oo /•oo 

= / / W{E,T)W{E' ,T)p{E, E')a{E)(7{E')dEdE' , (29) 
Jo Jo 

where —1 < p{E,E') < 1 is the energy dependent correlation coefficient. Notice that the variance 
depends on the correlation of our random function between two energies. If we naively propagated 
the uncertainty as the transform of the standard deviation: 

poo 

ax{T)= / W{E,T)a{E)dE, (30) 
Jo 

we would generally over-estimate the uncertainty, as seen in the quadrature difference between 
these two error assignments. 

/• oc /• oc 

al{T)-al{T)= j / W{E,T)W{E' ,T)[p{E,E') - l\(T{E)a{E')dEdE' (31) 
Jo Jo 

Since W{E,T) and (t{E) are positive definite and the quantity p{E,E') — 1 < 0, the difference is 
always less than or equal to zero. Thus, inclusion of these energy correlations reduces the total 
uncertainty in the thermal rates. What form these rates and errors take, depends on the what 
type of reaction we are dealing with and how we have assigned systematic uncertainties. Since 
we have treated the systematic errors as normalization errors independent of energy, it does not 
matter if we treat them in the integral or not. Actually, performing the integral both with and 
without including the systematic errors offers a nice way to double check the numerical integration. 
We now discuss the reactions important for primordial nucleosynthesis and their fits based on the 
former procedure. 
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III. RESULTS 



A. Cross Sections 

Keeping in mind our efforts to maintain a rigorous and model-independent analysis, we now 
implement this prescription for the set of nuclear reactions that are most important for big bang 
nucleosynthesis. Along with the neutron lifetime and Newton's Gat, eleven key nuclear reactions 
dominate the uncertainties in the BBN calculation of the light element abundances (Smith, Kawano 
and Malaney 1993 js^); determined by calculating the logarithmic derivative of the predicted 
abundances with respect to each of the reaction rates 37']. Thus, the choice of nuclear compilation 
with either its cross sections or thermal rates and their uncertainties, will determine the accuracy of 
the final predictions. The important work of Smith, Kawano and Malaney (1993) set a benchmark 
and their error budget has been the standard for Monte Carlo work. Nollett and Buries (2000) 
create their own compilation, but do not present portable fits of their cross sections and thermal 
rates. The NACRE collaboration (Angulo et ai; 1999 represents a large effort to critically 

evaluate the available nuclear data, presenting their adopted fits along with estimates of their 
uncertainties. Cyburt, Fields and Olive (2000) reanalyze a subset of the NACRE compilation in 
a simple, but uniform way in order to establish a more rigorous error assignment. Based on these 
most recent analyses, and the accuracy with which the WMAP satellite was able to determine the 
baryon density, it is clear that a rigorous and self-consistent prescription for dealing with nuclear 
data and deriving accurate representations and uncertainties must be established. It is with this 
main goal in mind, that we have developed the prescription in the previous section. 

There are two kinds of reactions, those induced by neutrons and those induced by charged 
particles. The cross sections for these reactions are generally decomposed into forms that behave 
more smoothly than the cross sections. In general, low energy cross sections scale with the square 
of the de Broglie wavelength, o" oc oc l/w^, where v is the relative velocity between the incident 
and target particles. There are further modifications to this behavior depending on the type 
of interactions involve. The neutron induced reactions feel only the strong nuclear force. The 
transmission probability of a neutron hitting this sharp potential surface is proportional to v, thus 
the neutron induced reactions can be written as follows: 

-(^) - S 

where R{E) is usually a smoothly varying function of center of mass energy, and constant 
at low energies. Na is Avogadro's number. The charge induced reactions feel the long range 
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electromagnetic force, with a transmission probability exponentially suppressed by the Sommerfeld 
parameter, (. The charge induced cross sections can be decomposed into 

a(E) = S(E)e.,(-2.0^ ^^^^ 

E 

where S{E) is the astrophysical S-factor, and the Sommerfeld parameter is defined by 

C.....4^J^V'^iftV'. (3.) 



2E J 271 \E 

Here the ZiS and /i are the charge numbers and reduced mass of the reactants, a is the fine- 
structure constant and Eg = 211"^ Z'fZ'^cP' is the Gamow energy. The S-factor, S{E), can also 
be a slowly varying function of energy. 

In the following we evaluate best-fits and uncertainties in R{E) and S{E), following our above 
statistical procedure. We use polynomial fitting functions, y{x) = J2n=o'^'^^^^ where the degree 
N of the polynomial is allowed vary until a minimum Xu is found. The data used in the following 
discussion has been gathered largely with the use of the NNDC's website 

1. n-decay and Newton's 

The lifetime of the neutron and Newton's G^r are key in determining the amount of ^He, being 
dependent on the neutron abundance at the deuterium bottleneck, they also dominate the ^He 
uncertainty. The lifetime of the neutron is key in determining the rate of neutron-proton inter- 
conversion. Reactions such asn+Ve p'^+e~ and n p'^ +e~ + i>e, have a common normalization 
and thus can be scaled with the mean neutron lifetime. The propagation of the neutron lifetime 

n 

uncertainty into the light element abundance predictions was first explored by Olive et al. 56] and 



in subsequent works 



33 



34 



35 




36l. l37l. l38l. l39l. l4Cll| . Newton's Gjv enters into the BBN calculation 



through the universal expansion physics. The effect of the gravitational constant's uncertainty has 
been previously examined by Scherrer and agrees well with this work's results. We adop t the 
recommended neutron lifetime and Newton's Gat from the Particle Data Group (2002) |2a] with 
Tn = 885.7 ± 0.8 sec and Gat = (6.673 ± 0.010) x 10^*^ (cgs). 

2. p{n,j)d 

Knowing the p{n, reaction is key in determining the end of the deuterium bottleneck and 
thus the onset of big bang nucleosynthesis. This radiative capture reaction is measured sparsely 



14 




FIG. 1: The reaction rate data for p{n,'y)d. The soHd Une represents the best fit, whilst the dashed the 
l-sigma error bars. The fit is an R-matrix calculation by Hale & Johnson (2003) [Sg. The data is shown 
with their respective l-sigma error bars. 

in the energy range of interest for BBN, .01 — 1.0 MeV. It is because of this lack of data that we 
must rely on a constrained R-matrix fit using elastic p — p, n — p scattering, and both unpolarized 
and polarized j — d photo-disintegration data, in addition to the sparse np-capture data of Nagai 
;59>] and Suzuki 1995 S^l- We adopt the R-matrix calculation of Hale and Johnson (2003) 
who have used the data discussed above to determine the np-capture cross section and its 
energy dependent uncertainties. This information was graciously provided by G. Hale upon private 
communication. This rate is now know to better than 2.5%, about a factor of 2 improvement over 
previous studies. We do not calculate a systematic error for this reaction. 

3. d(p,7)3He 

The d{p, 7)^He reaction is the first in a chain of reactions that rapidly burn deuterium after 
the deuterium bottleneck into ^He and eventually ^He. There are few data sets for this reaction in 
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FIG. 2: The reaction rate data for d{p, 7)^He. The sohd hne represents the best fit, whilst the dashed the 
1-sigma error bars. The data is shown with their respective 1-sigma error bars. The bottom panel shows 
the residual scattering in the data about our best fit, where our errors are set to ±1. The light curves are 
the NoUett and Buries 133 best fit and la errors. 



the BBN energy range. We consider the data sets of Bailey (1970) |6l|, Griffiths (1962,1963) 



Ma (1997) 



1, Schmid (1995,1996) 3 

and Casella (2002) j65l |. Some of these data sets warrant 
detailed consideration. The Casella data is the most recent measurement of this cross section and 
serves to anchor the low energy behavior of this reaction. This data has not been included in 
older analyses, only in this and two more recent BBN compilations by Cuoco et al 
al. 0i- It has been suggested that the 1963 Griffiths and 1970 Bailey experiments used incorrect 
stopping powers, and thus their low energy behavior is ~ 15 % too high 



n 



6g] and Coc et 



63 



6J]. Since the Casella 



data dominates the low energy behavior of the cross section, inclusion of the Bailey data does not 
affect this region of the cross section, thus we find no reason to omit it from our analysis. The 
1963 Griffiths data however, does not have a clear discussion of the normalization uncertainties, 
thus we exclude this data set from our analysis. The Schmid data sets suffer from poor energy 
resolution, with typical uncertainties in energy greater than 10%, which have not been included in 
their cross section errors. We thus exclude the Schmid data sets from the analysis. The inclusion 
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of the new Casella data, greatly increases the accuracy of the cross section, when compared to 
previous analyses. 

The discrepancy systematic error is S^isc = 0.0345, the intrinsic normalization error is Snorm = 
0.0528 and the total systematic error is 5 = 0.0631. 

4. d{d,nfRe 
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FIG. 3: Same as fig.Elbut for d{d,n)^lle. 

The d{d, n)^He reaction is the dominant deuterium sink during primordial nucleosynthesis. We 
consider the data sets of Brown (1990) 0|, Krauss (1987) 

Q, Ganeev (1958) Q, Arnold (1954) 
0, McNeill (1951) Q, Research Group (1985) Q, Preston (1954) 3,], Jarmie (1985) 75] and 
Schulte (197 2) [t^. Of these data sets, inconsistancies in the Ganeev data set found on the NNDC 
website j55|, create difficulties when trying to separate systematic errors from the total errors 
presented (e.g. unphysical statistical errors), thus we exclude this data set. The Arnold data exists 
only as a smoothed data set. This smoothing will artificially increase this data set's weight on the 
fit, thus we exclude this data set. The high energy Schulte data helps smoothly interpolate the gap 
between it and the low energy data. One may notice that the fitted curve falls below a majority of 
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the data, a seemingly bad "chi-by-eye" fit. The eye is misleading in this case. Since we are treating 
the correlations between data points explicitly, it is important to understand its impact. As we 
determined in the previous section, the statistical uncertainty plays a larger role than the total 
uncertainty of a particular data point. When a data set has very small statistical uncertainties, it 
gets more weight when determining the fit. This is exactly what we are seeing here. Though the 
Brown and Research Group data have small normalization errors, their statistical errors are large 
when compared to the statistical uncertainties in the Krauss data, where a < 1 %. Even though 
the Krauss data has larger normalization uncertainties, its statistical errors are significantly smaller 
than other data sets, and thus the Krauss data dominates the low energy behavior of the fit. It 
is interesting to note that if we turn off the correlations between data points and adopt the total 
uncertainties as the representative errors, we reproduce the mean value of the Nollett and Buries 
curve. 

The discrepancy systematic error is Sdisc = 0.0369, the intrinsic normalization error is Snorm = 
0.0400 and the total systematic error is (5 = 0.0544. 



5. d{d,p)t 



The d{d,p)t reaction is very similar to its mirror d(d, n)^He reaction, both in shape and mag - 
nitude. We consider the data sets of Krauss (1987) ^69], Brown (1990) Q, Preston (1954) 

0, Ganeev (1958) Q, 



Arnold (1954) 



Neill (1951) 



7l|, Davenport (1953) |77|1, Research Group (1985) [73], Ganeev (1958) 



Mc- 



72l |. and Gruebler (1981) [7a|. We exclude the data sets of Ganeev and Arnold for 
the same reasons as for the d{d, n)'^He reaction. We again see the statistical uncertainties in the 
Krauss data pulling the fit below the Research Group and Brown data. Again, this is entirely due 
our explicit treatment of the correlations in the data. If we turn off the correlations and adopt the 
total uncertainty as the representative uncertainty, we again reproduce the mean value curve of 
Nollett and Buries. 

The discrepancy systematic error is Sdisc = 0.0487, the intrinsic normalization error is Snorm = 
0.0560 and the total systematic error is (5 = 0.0742. 



6. ^ile{n,p)t 

The ^He(n,p)t reaction is responsible for the inter-conversion of mass 3 elements, maintaining 
an equilibrium relation between the two elements while this rate is fast when compared to the 
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FIG. 4: Same as fig. |3 but for d{d,p)t. 



Hubble exp ansion rate. We consider the data sets of Brune (1999) , Costello (1970) |80|, Coon 
(1950) sj> Gibbons (1959) 3> Macklin (1965) 83], Batchelor (1955) Q, Borzakov (1982) Q 
and Alfimenkov (1980) 3]. We exclude the Costello data because of poor energy resolution, the 
Coon and Macklin data because of little or no error information, and the Alfimenkov data because 
the reference was not available. There is a lot of data above 1 MeV for this reaction. In order 
to fit all of the data, we would need many fitting parameters. Since the energy range relevant for 
BBN is below 1 MeV, we do not use data above 1 MeV. As one can see, the fit is dominated by 
the Brune data. 

The discrepancy systematic error is 5 disc = 0.00703, the intrinsic normalization error is Snorm = 
0.0468 and the total systematic error is (5 = 0.0473. 



7. t{d, n)4He 



The t{d, n] 



Allan (1951) 



He reaction is a main production route to ^He. We will consider the data sets of 



Conner (1952) |2l|, Davidenko (1957) Q, Jarmie (1984) ^] and Research Group (1985) ^]. We 
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FIG. 5: Same as fig. 13 but for ^lie{n,p)t. 



exclude the Allan and Argo data sets because of uncertain normalization error assignments, and 
the Davidenko data set because the reference was not available. We again exclude the Arnold data 
set because of their smoothing their data. The Conner data assumes the cross section is isotropic. 
This assumption is good up to energies of ~ 240 keV, thus we exclude any Conner data that lie 
beyond this energy. 

The discrepancy systematic error is Sdisc — 0.0218, the intrinsic normalization error is dnorm — 
0.0401 and the total systematic error is 5 = 0.0456. 



8. 3He(d,p)4He 



The ^He(d, p)^He reaction is also a main route for producing ^He. We consider the data sets 
of Arnold (1954) [ZH], Bonner (1952) jaJLGeist (2000) Q, Krauss (1987) Q, Kunz (1955) Q, 
Moller (1980) and Zhichang (1977) [23]. We exclude the Arnold data again, because of their 
smoothing the data. 

The discrepancy systematic error is Sdisc — 0.0268, the intrinsic normalization error is dnorm — 
0.0605 and the total systematic error is 5 = 0.0662. 
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FIG. 6: Same as fig.Hbut for t{d,n)*}ie. 



9. 3He(a,7)^Be 



The ^He(a,7)''Be reaction is responsible for the production of ^Li in a high baryon density 
(?7 > 3x10"^°) universe. Its uncertainty dominates the prediction of the ^ Li abundance predi ction . 



We consider the data sets of Holmgren ( 1959 ) 



Krawinkel (1982) 



in2l |. Robertson (1983) ^M^, Hilgemeier (1988) 



9g|, Parker (1963 ) linnll . Nagatani (1969) 



1041 1 and Osborne (1984) 



ioi|, 



Following the suggestion of Hilgemeier, we renormalize the Krawinkel data by the factor 1.4, 
correcting the helium gas density. 

The discrepancy systematic error is 5 disc — 0.1482, the intrinsic normalization error is Snorm — 
0.0814 and the total systematic error is (5 = 0.1691. 

This reaction is also very important for stellar physics, in particular neutrino production. The 
low energy behavior of this reaction rate determines the flux of '^Be and neutrinos coming from 
the Sun. We believe it is inappropriate to base the low energy value on an average of extrapolated 
points, and recommend our adopted method of a global analysis of the data and its uncertainties 
and then extrapolating a low energy value. We get a value of 5*34(0) = (1.0ib0.169)(0.386±0.020) = 
0.386 zb 0.068 keV b for the astrophysical S-factor. This is significantly lower than the values 
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FIG. 7: Same as fig. H but for 3He(d,p)4He. 



determined by Adelberger et al. S^f^ = 0.53 ± 0.05 keV b jlOgl . the NACRE collaboration 
^.nacre ^ ^ q qq j^^y ^ Q ^^^^ ^^le Cyburt, Fields and Olive Q renormalized NACRE rate 
5*34^^ = 0.50 lb 0.05 keV b, though all determinations are based primarily on the same data. The 
Osborne data dominates the fit at low energy, causing the downward turn of the S'-factor. This 
turn is also seen in the Nollett and Buries compilation The model independent approach 

adopted in this work and in the work of Nollett and Buries should not be used for extrapolation, 
as these methods are meant to describe the data alone and thus are only valid where data exists. 
However, the inclusion of any theory fitting this data will still have to include the systematic errors 
similar to the ones discussed in this work. More measurements with E < 0.5 MeV will be able to 
more precisely determine 534(0). 



10. i(a,7)'^Li 



The i(a,7)^Li reaction is important for ^Li production in a low baryon density (77 < 3 x 10 ^'^) 
universe. Its uncertainty do minat es the theory prediction of ^Li's abundance here. We consider 
the data sets of Brune (1994) [mal, Burzynski (1987) JM], Griffiths (1961) [liia], Holmgren (1959) 
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FIG. 8: Same as fig. 13 but for 3He(a,7)^Be. 



, Schroder (1987) 



110| and Utsunomiya (1990) »llll. nomiya data set 

because of the lack of a normahzation error discussion. Smith, Kawano and Malaney and 
Nollett and Buries make the point that these Coulomb-breakup measurements are not yet 
reliable as this process is not yet completely understood, thus making the case for new experiments 
to be performed with E < 0.2 MeV. 

It is clear that the Holmgren and Schroder data are far from the best fit curve, outside of 
their assigned normalization errors. The visible discrepancy is forcing the sytematic error to be 
quite large. The Holmgren data also pulls the S'-factor fit down at ~ 0.6 MeV. If reason, other 
than the visible discrepancy exists to exclude these data, the fit would be dominated by the high 
precision Brune data with an overall 6% normalization error. 

The discrepancy systematic error is 6disc = 0.1788, the intrinsic normalization error is Snorm = 



0.1468 and the total systematic error is 5 = 0.2313. 
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FIG. 9: Same as fig. H but for i(a,7)^Li. 



11. '^Be{n,pyU 



The ''Be(n,p)^Li reaction is responsible for the inter-conversion of mass 7 elements at high 
baryon density [r] > 3 x 10 "^*^) This reaction has only one data set in the exoergic direction. The 
data set of Koehler (1988) |112l |. This data set does not extend very far into the energy range of 
interest for BBN. We must rely on the data for the endoergic re verse reaction, ^Li(p, n)^Be . W e 



consider the data sets of Gibbons (1959) Sekharan (1976) \UM and Taschek (1948) H- 
We use the principal of detailed balance to transfo rm t he '^Li(p, n)'^Be data into ''Be(n,p) ''Li data. 



Using t he Q- value from Audi and Wapstra (1995) 



website 



1151 ] available at the US Nuclear Data Program 



llfil ]. Q = 1.644168 ±0.000668 MeV. We ignore the lowest energy points derived from the 



reverse rate as they are sensitive to the precise value of Q, ignoring values that change significantly 
when Q is varied within its uncertainties. We should note that the Koehler data extends down 
to well below 1 eV, we choose not to show the data as its roughly constant and to emphasize the 
energy range important for primordial nucleosynthesis. 

The discrepancy systematic error is Sdisc = 0.0159, the intrinsic normalization error is Snorm = 
0.0448 and the total systematic error is (5 = 0.0475. 
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FIG. 10: Same as fig.Hbut for ^Bc(n,p)^Li. 



12. 



The ^Li(p, a)^He reaction is the dominant destruction chann el of ''Li at low baryon den sities 
(t/ < 3 X 10"^°). We consider the data sets of En gstler (1992) [113], Harmon (1989) [ud], Lee 
(1969) [ngl, Rolfs (1986) [ij and Spinka (1971) [12l|. We exclude the Harmon data because 



it is based on a measurement relative to ^Li(p, a)^He at energies E > 150 keV by Shinozuka et 
al (1979) |l22| . All but 3 points lie below this energy range, thus this measurement relative to 
^Li(p, a)^He is not valid at these energies. One may consider using the 3 points that are measured 
at appropriate energies, but it does not change our fit significantly. We also exclude the Lee data 
set as the reference was unavailable. 

This reaction has the largest Gamow energy of all the reactions we consider, and thus is the most 
susceptible to electron screening effects. In fact, the low energy behavior of this reaction is modified 
by electron-screening effects in the experimental set-up. This behavior can be parameterized as 
'7exp(-£') = c"barc(^ + Ue), relating the experimentally measured cross section to the bare nuclear 
cross section (i.e. no electron screening), where Ue is the screening potential ,123]. Where Ue is a 
measure of how much the Coulomb barrier has been reduced due to electrons screening the bare 
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FIG. 11: Same as fig. |21 but for ^Li(p, Q!)^He. Also shown is the e-screening corrected S'-factor and its 
uncertainty. 

nucleus (e.g. ''Li). For the experiments that screening is important (i.e. Engstler), Ue = 245 it 45 
eV is found to be the best fit. This agrees with the determinations of Engstler et al, who find 
Ue = 300ibl60 eV using an approximation of the e-screening correction to the observed cross section 
shown above. Englster used the high energy data to determine the best fit, and then extrapolate 
this to determine the screening potential. We have fit all the data, including the electron screening 
potential self-consistantly. 

The discrepancy systematic error is 6disc = 0.0194, the intrinsic normalization error is Snorm = 
0.0769 and the total systematic error is 5 = 0.0793. There is an additional energy dependent 
systematic error induced because of our electron screening correction. This can be well accounted 
for by: 

6e-scr = 0.02015 {^^) exp (-15.34^), (35) 
again added in quadrature with the other systematics. 
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B. Thermal Rates 



For the BBN temperature range, Maxweh-Boltzmann phase-space distributions are an excellent 
choice for baryons (eqn. 125(1 and the thermal rates become, 

where (y{E) is the cross section, not the standard deviation. For neutron induced reactions, using 
eqn. we find 

2 /""^ / E \ 



and for charge induced reactions, using eqn. 1331 we find 

■ 8 ^ 



E fE\^'^ 



kT \ E 



dE. (38) 



It is trivial to determine the weighting functions with these relations. For the neutron- and charge- 
induced reactions the respective weighting functions are: 



8 



1/2 



E fEg\^'^' 



kT \ E 



(40) 



After these integrals are performed numerically, we must find some representation of these 
thermal rates to implement into the BBN code. We will look at some cases in which the above 
integrals can be done analytically. This will ultimately guide us in determining the functional 
forms for these rates. 

Typically, neutron induced reactions follow the behavior noted in the previous chapter and 
are particularly smooth over the data energy range coverage. Thus a simple polynomial in E^^'^ 
will generally suffice. In this case, the integral can be performed analytically and the numerical 
integration serves as a test of the integrator. 

Most of the reactions are non-resonant charge induced reactions. In this case the integral 
cannot be done analytically, particularly at the temperature ranges relevant for BBN. In order to 
understand the reason for this and gain some insight for a possible functional form, we will look at 
the case where the temperature is much smaller than the Gamow energy, kT <^ Eg. We would like 
to turn the weighting function in eqn 001 into something more familiar, like a gaussian. To do this, 
we Taylor expand the argument of the exponent, about some energy Eq, defined such that the first 



27 



derivative with respect to energy is zero at Eq. Doing this we find, Eo/Eg = {kT/2Egf/^. Also 
needed for this analysis is the width of this gaussian, which we can find by evaluating the second 
derivative of the exponent argument at £"0. We find, do/Eg = 2/\/3{kT/2Eg)^/^ . The relevant 
perturbative parameter is the ratio of this width to effective energy, uq/Eq = {^2kT /21EgY/^ . 
In order for this gaussian approximation of the integrand to converge, the width must be much 
smaller than the effective energy. Not only does this demand kT <C Eg, but since the power is 
smah {kT/EgY/^ must be small. To get 10% convergence, kT < 10 ^Eg, which is well below the 
range relevant for BBN. However, if we continue, we find that the rate transforms into: 



27 Eg ^ 
AkT 



SeffiEo). (41) 



If S{E) is taken to be a polynomial in E, then S^ffiEo) is a polynomial in (o"o/i?o)^ kT^^'^, 
where only the even powers of our perturbative parameter appear due to the symmetry of the 
gaussian. We adopt this form, allowing the order of the polynomial describing S'e//(-E'o) to vary as 
needed until an accurate fit is reached. Reactions with broad resonances modify the above form, 
with the Breit-Wigner form, with E = Eq; 

Seff{Eo) 



AkT 



l + {{Eo-Ej,)/rR/2)^' 



(42) 



where E^ and T^ are the resonance parameters and Seff{EQ) is a polynomial in powers of kT^/^. 

Reactions with narrow resonances are typically the sum of a non-resonant piece and a Breit- 
Wigner form. Since the resonance is narrow, we can treat the Breit-Wigner form as a delta function. 
In the case of a neutron induced reaction, the resonant part of the rate becomes: 

Kes = ^^^R{Er) exp i-En/kT). (43) 

We wish to reiterate here, that the cross section fits and their energy dependent uncertainties 
are numerically integrated. These exact results are subsequently cast into a usable form, fit to one 
of the forms mentioned above to within 0.1%. 

In thermonuclear rate compilations, such as this one, it is important to realize that most com- 
pilations rely on the same experimental data to derive their representations. Because of this, they 
must generally agree with each other over the range of validity. This range is shown it table ^ 
The limits are solely based on the maximum energy of the data used {kTmax ~ Emax/'i)- 

We first compare to the BBN reaction rate standard. Smith, Kawano and Malaney (1993) [3^ . 
As we see from fig. I12( there is overall agreement with our compilation and theirs. The curves tend 
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TABLE II: This table shows to what maximum temperature each thermal rate is valid for. This is solely 
due to the data used in the analysis. Also shown are the systematic errors. 



Reactions 
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0.05oO 
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^lle{n,p)t 
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0.0071 


0.0468 


0.0473 


t{d, n)4He 


2.3 


0.0218 


0.0401 


0.0456 


3He(d,p)''He 


2.3 


0.0268 


0.0605 


0.0662 


3He(a,7)^Be 


7.8 


0.1482 


0.0814 


0.1691 




3.9 


0.1788 


0.1468 


0.2313 


^Be(n,p)^Li 


11.7 


0.0159 


0.0448 


0.0475 


^Li(p, a^Re 


3.9 


0.0194 


0.0769 


0.0793 



to diverge at high temperature, where there is no data pinning down the high energy behavior. 

The disagreement with the rate for 7)'^He is almost entirely due to the use and exclusion of 

different data sets. The disagreement at low temperatures for ^Be(n,p)^Li, is most likely due to 

their taking a minimum energy when integrating this rate {Emin = 1 keV). We cannot compare 

I I 

directly to Nollett and Buries (2000) 38] as they do not present thermal rates, but differences are 
attributable to the differences in cross section representations already discussed. The "wiggles" 
seen in n)^He, ^He(d, p)'^He and ^Be(n,p)''Li are due to the slightly different values adopted 
for the resonance parameters of each reaction. 

We now compare to the work of Cyburt, Fields and Olive (2001) [s^, which used renormalized 
NACRE rates and an estimate of the errors. Again, since these compilations are based on most 
of the same nuclear data the rates should be similar, as seen in fig. 1131 An interesting point is 
that for a majority of reactions, the Cyburt, Fields and Olive error budget underestimates the 
errors, compared to this work. The intrinsic normalization error we have included is typically 
as important as the discrepancy systematic error, which Cyburt, Fields and Olive assumed to 
dominate the error budget. Again the high temperature portion of the curves tend to diverge, as 
there is no data pinning down the high energy behavior. We point out that the d{p, 7)^He rate as it 
is systematically lower than our rate. This is entirely due to the NACRE collaboration's inclusion of 
data we have excluded. Also evident are the differences in the d{d,n)'^}le and d{d,p)t rates, where 
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FIG. 12: The thermal reaction rate residuals of the Smith, Kawano & Malaney (1993) [3^ compilation 
plotted against temperature in units of 10^ K. The solid line shows how their mean value compares to this 
compilations mean value. The dashed lines correspond with the la errors. 

our compilation falls below the fits adopted by Cyburt, Fields and Olive. The differences between 
this compilation and others for these rates, is that we have included the correlations between data 
points. It is for this reason, that the low energy data that has small statistical, but large total 
error is dominating the fits, pulling the low energy cross section down slightly. 

The overall agreement between different rate compilations is quite reassuring. The biggest 
advantage to our compilation is we have explicit treatments for dealing with correlated data, and 
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FIG. 13: The thermal reaction rate residuals of the Cyburt, Fields & Olive (2001) |39| reanalysis of NACRE 
[54! plotted against temperature in units of 10^ K. The solid line shows how their mean value compares to 
this compilations mean value. The dashed lines correspond with the Icr errors. 



estimating systematic errors. These systematic errors dominate over the statistical uncertainties in 
all cases, for the temperature range important for BBN, T ~ (0.5 — 1.2) x 10^ K as seen in fig. 1141 
Properly treating the correlation between fitting parameters when propagating the cross section 
fits into thermal rates has a noticeable reduction in the statistical uncertainties. This reduction is 
maximized when there is a cross over between terms in the fit polynomial, when one term goes from 
being dominant to being sub-dominant and vice versa. It is also reassuring to see the statistical 
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uncertainty become the dominant contribution to the error at high temperature or energies, where 
there is no data. 
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FIG. 14: The thermal reaction rate relative errors for the 10 reactions fitted in this compilation. Tlie 
dashed and dotted curves are the thermally averaged statistical errors with and without treating the energy 
correlations. The dashed-dotted curve shows the total systematic errors, and the solid curve shows the total 
thermal error. 



32 



C. Light Element Predictions 



Adopting the thermonuclear reaction rates discussed in the previous section, we discuss their 
impact on BBN predictions and on the general concordance of the BBN predictions with the light 
element observations and the CMB. Shown in figureElare the light element abundance predictions 
and their percent errors using this work's nuclear compilation. 
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FIG. 15: Shown in the left panel are the light element predictions using this works nuclear rate compilation 
and uncertainties. The mass fraction of "^He (Vp) and the mole fractions, D/H, ■^Hc/H and ^Li/H are plotted 
against the baryon-to-photon ratio. The width of each curve represents the 1 a errors in the light element 
predictions. The right panel shows the relative uncertainties in percent of the light element predictions. 



Before we examine the concordance between this compilation's predictions and observations, 
we should verify the agreement of previous compilations and qualify their differences. As has been 
discussed in previous chapters, the BBN compilations of Smith, Kawano and Malaney |3J], Nollett 



and Buries [3g| and Cyburt, Fields and Olive 



39|, should all roughly agree as they are largely 



based on the same nuclear data. Any differences in their predictions will arise entirely from each 
compilation's derivation of reaction rates and their uncertainties, and the data each uses. 

Since Cyburt, Fields and Olive and Nollett and Buries both compare directly to Smith, Kawano 
and Malaney showing rough agreement, we choose to compare only to the former two compilations. 
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Plotted in figure is the residual between Cyburt, Fields and Olive |39| and this compilation, 
where zero and ±1 represent the means and standard deviations of this compilation. The '^He 
mass fraction (Yp) is in good agreement with this compilation. The ^He abundance error is slightly 
increased due to the inclusion of the uncertainty in Newtion's Gat. This compilations treatment of 
the data, leads to differences for the D, '^He and ^Li yields. On the high t] side {rj > 3 x 10^^'^) the 
changes are due to the d{p, 7)^He reaction. The new cross section is larger than the one determined 
by the NACRE-based compilation jsJ] of Cyburt, Fields and Olive, causing a subsequent drop in D 
yields with a simultaneous jump in ^He and ^Li yields. For low values of rj, only the mean value of 
^Li is significantly different. This is due entirely to a slightly lower cross section for the t(a,7)^Li 
reaction used here. The errors of the Cyburt, Fields and Olive compilation are generally smaller 
than this compilation's errors. This is due to the fact that this compilation has an additional 
intrinsic normalization error added in quadrature with the discrepancy normalization error. 




FIG. 16: This figure shows the difference between light element yields using this compilation and that of 
Cyburt, Fields and Olive [s^. The solid and short-dashed curves show the Cyburt, Fields and Olive yield 
means and standard deviations with respect to this compilation's means and standard deviations, seen here 
as zero and ±1, respectively. 

n 

Similarly, plotted in figure El is the residual between Nollett and Buries |3^] and this compi- 
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lation. The central values of the abundance yields are in good agreement. The most noticeable 
difference is in the ''Li yields, with this compilation having slightly lower values. The high rj differ- 
ence is due to differences in the reactions d{p,j)^ILe and ^He(d,p)^He, while on the low rj side the 
differences are primarily due to the reaction t{d, n)^He. The errors of Nollett and Burle's compila- 
tion are comparable, but generally smaller than this compilation's errors, because of the statistical 
nature of their errors. The larger ^He error is due to the larger uncertainty in the neutron lifetime 
adopted by Nollett and Buries. 




FIG. 17: This figure shows the difference between light element yields using this compilation and that of 
Nollett and Buries js.^- The solid and short-dashed curves show the Nollett and Buries yield means and 
standard deviations with respect to this compilation's means and standard deviations, seen here as zero and 
±1, respectively. 



It is reassuring that this new BEN nuclear compilation a gree s quite well with the previous studies 
of Smith, Kawano and Malaney |34i]. Nollett and Buries jSa ] and Cyburt, Fields and Olive 39 1, 
as well as the two more recent calculations by Cuoco et al. and Coc et al. 67], though the 
new compilations do not present rate representations that can be compared directly. This work's 
rigorous treatment of systematic uncertainties also suggests that these previous works may have 
underestimated the true error budget in the light element abundance predictions. This comparison 
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37|. They can be 



also suggests where new nuclear data will be most useful. To illustrate this point, scalings are 
created that explicitly show how the baryon-to-photon ratio rj, Newton's G and the reaction rates 
affect the light element abundance predictions. These scalings are calculated numerically, by finding 
the logarithmic derivatives of the light element abundance predictions with res pec t to parameters 
and key reaction rates, relative to a fiducial model where rj = 6.14 x 10^^'' 
used to either predict light element abundances or propagate uncertainties, but these scalings are 
only approximate and will change for models with rj very far from its fiducial value. We use them 
here, only to discuss how BBN's predictions depend on the various inputs. The nuclear reactions 
are parametrized here through Ri, where i refers to the subsection number assignment for that 
reaction in sect. IIIlXl fi.e. R2, R4, and correspond with the p(n,^)d, d{d,n)'^Iie, and d{d,p)t 
reactions respectively). The Ri can be thought of as reaction normalizations, such that the current 
compilation is Ri = 1.0. The scalings are: 

/inlO \0-39 / \ 0.72 / p \0.35 

Yp = 0.24849 ( ^ ) ( ^ ) (^] i?ooo6i?o.oo5^o.oo5 (44) 

V 6.14 J \Tn.oJ \Gn,oJ 



n /TfllO \^l-62/ \0.41/ pi \0.95 

'»S-^-^^Kw) (^) {§^) l^i'-W^Rf-- (45) 

3tt /inlO \-0.59/ \0.15/ p \0.34 

^e_tie ^ 10 086 ^ \~) \ d-o.77 dO.ss D-0.25 0-20 D-0.17 dO.os 



10^^ = 10.086 ^f-f — V''i?^-^«i?5-'^-"^i?4-^X ^2 (46) 



7t ■ /inlO \212 / X 0.44 / p \-0.72 

IQioJ^ = 4.364(^) (^) (^] i?i-3X''V''V-''<''«3''V''-(47) 
H \ 6.14 J \Tn,oJ \Gn,o J 

As clearly seen in the scalings, Yp is dominated by the neutron mean lifetime and Newton's 
Gat, while the reactions n)^He, d(d,p)t and p{n,^)d only slightly change the predictions. 
The dramatic drop in sensitivity (scaling powers ~ 0.5 to ~ 0.005) is seen in the other light 
element scalings, meaning those reactions do not contribute significantly to the overall theory 
predicions. Thus, for brevity these reactions are left out of the scalings for D, ^He and "^Li. 
For an accurate D prediction, d{d,n)^He and d{d,p)t are key, with d{p,^)^ile following close 
behind. ^He is the least sensitive to which nuclear compilation is used, though improvements 
in its prediction propagates into an improved '^Li prediction. For high baryon densities (?? > 
3 X 10^^^), the reactions ^He(a,7)^Be and ^He((i,p)^He dominate the ^Li predictions, whilst for 
low baryon densities [r] < 3x lO"^'') their mirror reactions are dominant, t(a,7)^Li and n)^He. 
With the precision of p{n,j)d being < 2.5%, it only enters at the percent or sub-percent level 
in the light element prediction uncertainties, and thus is not the dominant error. With this new 
nuclear compilation and its error budget we are well poised to test the overall concordance between 
primordial nucleosynthesis' predictions, the observations of the light element abundances and of 
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the CMB anisotropy. 



IV. DISCUSSION 



A. Light Element Observations 



With the hght element predictions of D, ^He, ^He and ^Li in hand, we set out to compare them 
directly to observations. Deuterium is measured in high-redshift QSO absorption line systems 
via its isotopic shift from hydrogen. Under the wel l-fou nded assumption that the only significant 
astrophysical source of deuterium is the big bang 
depletion in these high-shift systems to be less than 1%. Thus, making D nearly primordial and a 
direct probe of big bang nucleosynthesis. In several absorbers of moderate column density (Lyman- 
limit systems), D has been obs erved in multiple Lyman transitions. We adopt the two deuterium 



12J], one can estimate that the amount of D 



values from Kirkman et al. 



127l |. one being the world av erage of the 5 



3est deuterium measurements 



including both single and multiple absorption systems 



125, 



12e, 



12l 



12H : 



(2.781°:^^) X 10-'. 



(48) 



and second, the average of the 2 multiple absorption line systems 



126, 



127|; 



(2.49t°:20) X 10"^ 



(49) 



As noted in Kirkman et ai, the P^r degree of freedom is rather poor for the world average 
D value {x1 = 4.1). Many possibilities exist that can ex plain this poor x^, underestimated errors, 



129l |. However, since we are only dealing 



correlations with column density, and other systematics 
with 5 systems, any of these conclusions can be reached. The 2 multiple absorption line systems 
agree quite well with each other, however this could also be due to low number statistics. Future 
observations will help address these concerns. 

Unlike D, ^He is made in stars, and thus co-produced with heavy elements. Hence the best 
sites for determining the primordial ^He abundance are in metal-poor regions of hot, ionized gas 
in nearby external galaxies (extragalactic HII regions). Helium indeed shows a linear correlation 
with metallicity in these systems, and the extrapolation to z e ro metall icitv gives the primordial 
abundance (baryonic mass fraction, Yp = pi^^/ps) |l3ol . Il3l1 . 1132I . [l3^ . We cite the 2 following 



values, as have 



3£, 



46 



3 



Yp = 0.238 ± 0.002 ± 0.005. 



(50) 
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Yp = 0.244 ± 0.002 ± 0.005. (51) 

determined by Fields and Olive (1998) and Izotov and Thuan (1998) 'l3lj | , respectively from a 

large body of data representing dozens of extragalactic HII regions. The difference between the two 
values is due primarily to adopted analysis techniques treating ^He emission lines and underlying 
stellar absorption, as most of the systems are the same between the two sets. Here, the first error 
is statistical and reflects the large sample of systems, whilst the second error is systematic and 
dominates. Since Izotov and Thuan do not quantify a systematic error, we a dopt the systematic 
error discussed in Fields and Olive and explored in Olive and Skillman (2001) 
by Olive and Skillman, these systematics need to be further explored. 

Helium-3 is observed as well, but through its hyperfine emission in the radio band, limiting 
observations to Galactic HII regions. The sam ple si ze of the ^He data is rather sparce and localized 
around a fairly narrow band in metallicities 



13J]. As suggested 



140|. Combined with a considerable dispersion, a 
model independent determination of the primordial ^He abundance is prohibitive. The Galactic 
evolution of ^He is also poorly understood, as it is not known if ^He increases or decreases from 
its primordial value |l4l| . manefesting itself as a large extrapolation error in model-dependent 
approaches. We thus, do not use ^He observations to probe primordial nucleosynthesis. 

The primordial '^Li abundance is determined from observations of old metal-poor stars, particu- 
larly those in the Galactic stellar halo (Population II). For very low metallicities, the ^Li abundance 



1421 ■ F^om this, a primordial abun- 



is found to be nearly constant, the so-called "Spite Plateau" 
dance is inferred . An analysis of a set of Pop. II stars with high signal to noise data was performed 
by Ryan et al. 143l |. taking into account various chemical and stellar evolution effects. Their 



primordial Li abundance is: 

'7 



—\ = (1.23 ± 0.06l[!-|) X 10-10 (95% CL), (52) 
where the small statistic al er ror is overshadowed by systematic uncertainties. A recent determina- 



tion by Bonifacio et al. jl44l |. based on observations of stars in a globular cluster, yields slightly 
different results. The difference is mainly attributable to the different methods used to calli- 
brate stellar atmosphere parameters, in particular the effective temperature. Their analysis yields 
^Li/Hp = (2.19lQ'3g) ^ IQ-i'^. The difference between these numbers is a measure of the systematic 
error, which has apparently been underestimated by Ryan et al.. We thus adopt both observations 
for use as probes of primordial nucleosynthesis. 
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Since standard primordial nucleosynthesis is a one parameter theory, depending on the baryon- 
to-photon ratio, rj or equivalently the baryon density, we can use light element abundance deter- 
minations to measure the baryon content of the universe. We discuss the implications of adopting 
the observations mentioned, on BBN concordance. 
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FIG. 18: Shown in the figure are the light element predictions. The mass fraction of ^He (Yp) and the 
mole fractions relative to hydrogen, D/H, '^He/H and ^Li/H are plotted against the baryon-to-photon ratio. 
The width of each curve represents the 1 cr or 68% confidence errors in the light element predictions. The 
outlined boxes represent the light element observational constraints on the baryon density. 



Shown in fig. are the light element predictions with outlined boxes showing the observational 
constraints and the r] ranges allowed by each. There is no value of baryon density for which any 
three abundance observations agree well, as seen quantitatively in tab. lIIII Treating all observations 
equally, we can only reliably constrain the baryon-to-photon ratio to lie between 1 < W^^t] < 7. 
There is only marginal agreement at the 95% confidence level. This marginal concordance is also 
eva l uated in previ ous works js^ . |3] with the use of theory and observationally based likelihoods 



56 



145 



Me, 



143]. If we limit ourselves to D only constraints we find that W^^r] = 6.28toit 

and 



5.92l;Q5g, for the multiple absorption and world averages respectively. 

This tension could either be pointing out unknown systematics in the abundance observations, 
or be telling us that there is new physics to be learnt. We address both of these in the following 
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TABLE III: This table lists the baryon density constraints placed by various light element observations 
using this works theory predictions. For comparison, the WMAP team's result for the baryon density is also 
shown. The numbers cited are the mostly likely values and their respective 68% central confidence limits. 
Since the 144] ^Li constraint lies above the "dip" in the theory prediction, it has two distinct predictions 
for the baryon density, a low baryon density constraint / and a high baryon density constraint //. 



Observations 
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0.0229 ±0.0013 
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0.0216l°:°g^? 


Yp = 0.238 ± 0.002 ± 0.005 [132] 






Yp = 0.244 ± 0.002 ± 0.005 fl33] 




0144+" °i2^ 

U.Ui44:_g QQgQ 
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U.UiiD_g QQ44 
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"4.4it°:^I 


u.uuoo_o 0008 

U.UiDi_o.o019 


WMAP (2003) iJJj 


6.14 ± 0.25 


0.0224 ±0.0009 



sections. An independent measure of the baryon density will eliminate it as a free parameter for 
BBN. This independent determination will act as a tie-breaker among light element observations 
and lead the way to understanding this tension, whether the disagreement results from underlying 
systematics or new physics. 

B. Observational Concordance 

As mentioned at the beginning of the work, the CMB anisotropies detail information about the 
shape and content of our universe. With the first data release of the WMAP team pj], several 
cosmological parameters have been measured to unprecedented accuracy, including the baryon 
density, which is measured to be, 

nsh^ = 0.0224 ± 0.0009. (53) 

This corresponds with a baryon-to-photon ratio oi r] = (6.14 ± 0.25) x 10^^'^. This is a 4% 
measurement, which makes it a sharper baryon probe than any light element currently is. Since 
we no longer are required to use the light element abundances to tell us the baryon content of 
the universe, the analysis completely changes. Now we can predict the light element abundances, 
with this baryon density and compare those predictions with the light element observations. With 
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WMAP's bar yon density we get: 



D/H 
^He/H 
^Li/H 



0.2485 ± 0.0005 



(2.55t°|J) X 10-^ 



;i0.12lO|7) X 10- 



(4.26+°!^) X 10-^° 



(54) 
(55) 
(56) 
(57) 



for the light element predictions with our new nuclear reaction network. Figure shows the 
predictions and compares them directly with observations. 
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FIG. 19: Primordial light element abundances as predicted by BBN and WMAP (dark shaded regions). 



region shows D/H = (2.78l^;||) x 10 



10- 



126, 



125 


12fi 


127 


128 



1231, while the dashed curve shows D/H = {2.AQZ^) ii) ^ 



127 1; (b) no observations plotted; see text (c) the light sh aded region shows Yp — 0.238 ± 0.002 ± 



0.005 !l32|, while the dashed curve shows = 0.244 ±0.002 ±0.005 Il33l; (d) the light shaded region shows 



l~l 



^Li/H = (1.23;° jg) X 10-1° |143|, while the dashed curve shows ^Li/H = (2.19;°;3^) x 10"!° |l44 



In order to quantify the level of concordance, we define an effective x ■ 



^2 



obs "wmap 
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where Aobs and A^^ap are the most Hkely values of the Hght element abundances for the adopted 
observations and that predicted with the WMAP baryon density, aobs and (Twmap are the corre- 
sponding 68% confidence errors. The values are shown in tab. IIVI 

TABLE IV: This table lists the effective x^'s for each observational constraint of the light element abun- 
dances, given the WMAP baryon density and this compilation's BBN theory. A Xe^y value smaller than 
unity means concordance, while a value large than unity shows discordance. The magnitude of discordance 
is measured by yXe//' ^ measure of how many "cr" of discordance exists. 



Observations 






D/H = (2.491^:2^) X 10-5 [126. 1271 


0.045 


0.212 


D/H = (2.78l[Jj^) X 10-5 [n^V2^\22^\M 


0.281 


0.530 


Yp = 0.238 ± 0.002 ± 0.005 [132] 


3.77 


1.94 


Yp = 0.244 ± 0.002 ± 0.005 [ml 


0.692 


0.832 


Yp = 0.238 ± 0.002 \132] 


25.9 


5.09 


Yp = 0.244 ± 0.002 [133] 


4.77 


2.18 


^Li/H = {1.23±Omtoit) X in-^°f143l 


10.72 


3.27 


^Li/H = (2.19t°:;^^) X 10-i°[144] 


4.50 


2.12 



As one can see, the two adopted observational values of deuterium agree with the BBN+WMAP 
prediction, both having Xe//'s smaller than unity. It is unclear if the slightly worse of the world 
average is due to unknown systematics or just poor statistics. Hopefully, with future automated 
searches, many more of these special absorption systems can be found. It is interesting to note 
that the WMAP baryon density contributes significantly to the uncertainty in the predicted D 
abundance. Future CMB experiments will reduce this uncertainty, at which time the BBN nuclear 
uncertainties will totally dominate the theory predictions. Thus motivating renewed efforts for new 
cross section measurements. 

The Izotov and Thuan value, Yp = 0.244 133l | value agrees with theory predictions only if the 



systematic errors are taken into account as discussed earlier. If they are ignored, this nu mber 
shows discordance at more than the 2-a level. The Fields and Olive value, Yp = 0.238 
shows discordance at the 2-a level with systematic uncertainties. If they are ignored here, the 
discordance becomes a 5-cj deviation. It is clear that a more detailed study of these systematics, 
including the effects of underlying stellar absorption and varying treatments of emission lines, 
is needed 134l |. One may also consider the new evaluation by Izotov and Thuan 135l |. finding 
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Yp = 0.2421 ± 0.0021, a nearly 3-cj deviation with this compilation's CMB+BBN predictions. We 
beUeve a proper accounting of the systematic errors will alleviate this discordance. 

The CMB itself, is also sensiti ve to the value of Yp. First attempts at co nstr aining Yp have been 



performed by Trotta and Hansen 



136l | and by Huey, Cyburt and Wandelt 



137i | . Future parameter 



studies, should include 1^ as a free parameter, rather than adopting the canonical value of 0.24. 
The CMB constraint offers an independent determination, free of the systematics plaquing the 
determination from extra-galactic HII regions. It also is a direct probe of Yp, so no extrapolations 
to zero metallicity are needed for the determination, just high precision CMB anisotropy data. With 
this data, we can use the CMB-determined Yp to quantify the level of observational systematics 
discussed above. By comb ining BBN predictions with CMB observations, we can also learn about 



stellar evolution Il36 



139|. One should also be mindful of the baryon density dependence on Yp 



given in eqn. Q especially when combining BBN and CMB results. 

The WMAP+BBN prediction for ^Li disa grees with bot h obs ervationally-based primordial ''Li 
abundances, with the Ryan et al. and Bonifacio et al. 144 1 numbers showing discordance at 
the 3 and 2-a level. As already mentioned, the difference between these two sets of observations is 
a measure of the systematic error due to the different methods used. This is not large enough to 
account for all of the discrepency between the observation-based and predicted values. An often 
discussed possibility is the depletion of atmospheric ''Li. This possibility faces the strong constraint 
that the observed lithium abundances show extremely little dispersion, making it unlikely that 
stellar processes which depend on the temperature, mass, and rotation velocity of the star all 
destroy ^Li by t he s ame amount. Uniform depletion factors of order 0.2 dex (a factor of 1.6) have 



been discussed 



148 1 . It is clear that either (or both) the base-line abundances of ^Li have been 



poorly derived or stellar depletion is far more important than previously thought. Of course, it is 
possible that if systematic errors can be ruled out, a persistent discrepancy in ''Li could point to 
new physics. 



C. Implications for non-standard BBN 

With the goal of maintaining concordance with observations, we examine how sharply we can 
deviate from the standard model. Often the effect of new physics can be parameterized in terms of 
additional relativistic degrees of freedom during the epoch of primordial nucleosynthesis, usually 
expressed in terms of the effective number of neutrino species, Nj^^eff- Traditionally, D or ^Li 
observations were used to fix the baryon density and the ^He mass fraction was used to fix Ny^f,ff. 
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These limits are thoroughly described elsewhere jSTj, l4ll . l42l . |4g]. Moreover, as we have noted, 
the observed ^He appears lower than the WMAP+BBN value. This discrepancy is likely due to 
systematic errors, but could point to new physics. Until this situation is better understood, caution 
is in order. Fortunately, in the post-WMAP era, we can now use the CMB-determined baryon 
density fean. to remove it as a free param eter f r om B EN theory and use any or all abundance 



observations to constrain A^'j^^e// |4y, |4 



66 



149 



150|. In particular, we have computed the 



likelihood distributions for N^^^ff using the WMAP rj and several of the light element observations. 



0.02 0.03 




10-10 10-9 

baryon-to-photon ratio r/ 



FIG. 20: The light element predictions plotted against the baryon-to-photon ratio for different values for 
N^^eff- The light shaded region corresponds with N^^^ff = 2.0, the medium shaded with N^^^ff = 4.0 and 
the dark shaded with N^,eff — 3.0. 

To first gauge what elements are sensitive to Ni^^^ff, we have plotted the primordial abundance 
predictions for the standard case, Ny^(.ff = 3, and two non-standard cases, Ny^(.ff = 2,4 in fig. QUI 
As readily apparent, ^He is the most sensitive element. If we understood the underlying systematics 
with the ^He observations better, this would be the ideal choice for picking an observation to make 
the constraint. However, since we are unsure about ^He, we must move to another observation. We 
see that D is the next most sensitive, and most notably, with this new compilation, the differences 
between the 3 A^^^e/j's are clearly resolvable. With current theory and observation uncertainties, 
''Li is not very sensitive to the relativistic degrees of freedom at the high baryon densities that the 
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CMB prefers, thus making it not suitable for this analysis. 

To gauge the kinds of constraints we can place, we calculate both ^He and D constraints seen 
in tab. EI The ^He observations that appear systematically low in the standard case, N,y,eff = 3, 
pull down the most likely value of Nj^^eff to lie between 2 and 3 depending on which abundance 
we favor. With our concerns of systematics in ^He observations, we do not put much weight in 
the ^He-based constraints, but we do note that the Izotov and Thuan 133 1 Yp determination is 



in fair agreement with the Standard Model of Particle phyics value of N^^eff = 3.0, as long as 
systematic errors are taken into account. However, the N^^^fj = 3.0 D prediction is in accord 
with the CMB baryon density, thus yielding less pull away from N^^^ff = 3.0 in the non-standard 
model. It is very interesting that each observation's most likely Nj^^^ff lie on op posite sid es of 
the standard BEN value of 3. Using the multiple absorption line system average 126l . Il27l |. we 
find N^^eff = 2.78t[];^^. With the world average, we find N^^eff = 3.65l:\il- Even though D's 
dependence on Ny^f.ff is smaller than its dependence on the baryon density, a sufficiently accurate 
measurement (e.g. WMAP) of ry, will help make D a more accurate probe. Thus demanding both 
improved nuclear data and more D observations. 



TABLE V: This table lists the N^ ^ff constraints placed by various light element observations using this 
works theory predictions and the WMAP team's baryon density, JIb/i^ = 0.0224± 0.0009 \v\. 
cited are the mostly likely values and their respective 68% central confidence limits. 



The numbers 



Observations 




D/H = ('2.49+^ i2l X 10-5 [126. 1271 


7Q+0.87 


D/H = (2.781^:38) X 10"''^ f125. 126. 127. 128] 


■^•"3-1. 30 


Yp = 0.238 ± 0.002 ± 0.005 UM 


2 26+°-3'^ 


Yp = 0.244 ± 0.002 ± 0.005 [133] 
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V. CONCLUSIONS 

Primordial nucleosynthesis has entered a new era. With the precision observations of WMAP, 
the CMB has become the premier cosmic baryometer. The independent BBN and CMB predictions 
for rj are in good agreement (particularly when D is used in BBN), indicating that cosmology has 
passed a fundamental test. Moreover, this agreement allows us to use BBN in a new way, as the 
CMB removes 77 as a free parameter. One can then adopt the standard BBN predictions, and use 
^CMB to infer primordial abundances; by comparing these to light element abundances in different 

45 



settings, one gains new insight into the astrophysics of stars, HII regions, cosmic rays, and chemical 
evolution, to name a few examples. Alternately, WMAP transforms BBN into a sharper probe of 
new physics in the early universe; with r/cMB fixed, all of the light elements constrain non-standard 
nucleosynthesis, with N^^eS being one example. 

As BBN assumes a new role, much work remains to be done. To leverage the power of the 
WMAP precision requires the highest possible precision in light clement observations. Further 
improvements in the primordial D abundance can open the door to D as a powerful probe of 
early universe physics. Improved ^He observations can offer new insight into stellar and chemical 
evolution in the Galaxy. And perhaps most pressing, the WMAP prediction for primordial ^He 
and particularly ^Li are higher than the current observed abundances; it remains to be resolved 
what systematic effects (or new physics!) has led to this discrepancy. 

This work has always been motivated by the idea of precision cosmology. We have laid out a 
rigorous procedure for determining best fit parameters and their uncertainties. We have explicitly 
taken into account the correlations among data points and their normalization errors. We found it 
necessary to define two systematic uncertainties, one is a calculation of the inherent normalization 
of the data. The second is a measure of how well different data sets agree with each other. This 
work generally agrees with previous studies, except in some special cases as discussed. 

Using these updated nuclear inputs, we compute the new BBN abundance predictions, and 
quantitatively examine their concordance with observations. BBN theory uncertainties are domi- 
nated by the following reactions: n)^He, d{d,p)t, 7)^He, ^He(a, 7)'''Be and ^He((i, p)^He. 
Reducing BBN's uncertainties will allow stronger statements about concordance. Depending on 
what deuterium observations are adopted, one gets the following constraints on the baryon density: 
^Bh"^ = 0.0229±0.0013 or J^b/j^ = 0.0216l^;gg^? at 68% confidence. If we instead adopt the WMAP 
baryon density, we find the following constraints on the effective number of neutrinios during BBN: 
K,eff = 2.781°: yg or Nj^^gff — 3.65lj;;3Q at 68% confidence. Concerns over systematics in helium 
and lithium observations limit the confidence of the constraints derived from this data. Further 
exploration of these systematics, given new observational techniques and more detailed models, 
will be most beneficial in understanding and ultimately reducing their effects. Deuterium suffers 
from a small sample size; a larger sample size will not only improve statistics but also allow the 
examination of possible systematics. With new nuclear cross section data, light element abundance 
observations and the ever increasing resolution of the CMB anisotropy, tighter constraints can be 
placed on nuclear and particle astrophysics. 

In closing, it is impressive that our now-exquisite understanding of the universe at z ~ 1000 also 
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confirms our understanding of the universe at z ~ 10^^. This agreement lends great confidence in 
the soundness of the hot big bang cosmology, and impels our search deeper into the early universe. 
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APPENDIX A: CROSS SECTION FITS 

1. d(p,7)3He 

S{E) = 0.2268(1 + 22.05-E + 30.77£;2 - 9.919£;3) eV b, 

2. d(d,n)3He 

S{E) = 0.05067(1 + 7.5ME - A.22bE'^ + l.SOS-E^ - 0.2041E;'') MeV b 

3. d{d,p)t 

S{E) = 0.05115(1 + 4.685^ - 1.021^^) ^ 

4. 3He(n,p)t 

R{E) = 6.846 X 10^(1. - 0.46474331 1577^; - OlE'^^'^ - 0.206566636058^; + 02E 
+ 145.303829979£;3/2 _ 5i7.845305322£;2 _^ 1061.59032882^5/2 - 1232.39931680£^3 
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+ 748.452414743£;^/2 _ 184.417975062£;^) cm^ g~^s-^ 



5. t{d, n) 



S{E) = 24.19(1 + 3.453E - 40.16^;^ + 285.6E^ - 596AE^ + 407.1£;5)/(1 + 



where 



Er = 0.0482 MeV and Tr = 0.0806 MeV. 



6. 



SiE) 




where Er = 0.183 MeV and Tr = 0.256 MeV. 

7. 3He(a,7)^Be 

S{E) = 0.3861(1 + 0.8195^; - 2.194^2 + 1.419£;3 - 0.2780^;^) keV b 

8. t{a,jy Li 

S{E) = 0.08656(1 + 0.6442E - 7.597^;^ + 12.16^;^ - 5.336^;^) keV b 

9. ^Be(n,p)^Li 

R{E) = 4.7893 x 10^^(1 - 4.12682044152£;V2 + 3.10200988738^; 

+ 15.8164551655£;3/2 - 45.5822669937£;2 + 54.7133921087£;5/2 _ 34.7483784037£;3 

+ 11.3599443403£;^/2 _ i.49669812741E^) + 1.0553 x 10V(1. + {{E - ER^i)/{0.5rR^i)f) 

+ 2.0364 X 10V(1. + {{E - ER^2)/{0.5TR^2)f) cm^ g-^s"! 

where Er^i = 0.32 MeV, Fr^i = 0.20 MeV, Er^2 = 2.7 MeV, and Tr^2 = 1-9 MeV. 

10. ^Li(p, a)^He 

S{E) = 0.06068(1 + 3.174^; - 7.586^2 + 8.539^3 - 3.216E'') MeV b 
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APPENDIX B: THERMAL RATES 

In adopting these rates, it is a good idea to assume the rate constant at temperatures above 
which these fits are no longer vaHd (see table ^ . This prevents the artificial divergence of the 
abundances. 

1. p(n,7)d 

Na{(tv) = 4.40654e4*(l.+.0457518*t912-2.47101*t9+4.17185*t932 

-3.44553*t9*t9+1.72766*t9*t932-.546196*t9**3 

+.106066*t912*t9**3-.0115306*t9**4+.536436e-3*t912*t9**4) 

2. d(p,7)3He 

NA{crv) = 7.30909e+3*t9m23*ex(-3.7209/t913) 

*(l.-10.3497*t913+63.4315*t923-209.780*t9 

+432.557*t943-571.937*t953+497.303*t9*t9 

-284.936*t943*t9+106.863*t953*t9-25.7496*t9**3 

+3.81387*t913*t9**3-.313823*t923*t9**3+.0108908*t9**4) 

3. d{d,nfRe 

NA{crv) = 1.00749e+9*t9m23*ex(-4.2586/t913) 

*(1. -9.59015*t913+65.2448*t923-247.756*t9 

+596.231*t943-941.064*t953+980.076*t9*t9 

-643.032*t9*t943+211.982*t9*t953+29.0491*t9**3 

-66.1847*t913*t9**3+31.6452*t923*t9**3-7.15147*t9**4 

+.372749*t913*t9**4+.208645*t923*t9**4-.0545129*t9**5 

+.00536216*t913*t9**5-.000157984*t923*t9**5 

-.457514e-5*t9**6+2.123592e-9*t913*t9**6) 
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4. d{d,p)t 



Na{(Jv) = 3.91889e+8*t9m23*ex(-4.2586/t913) 
*(l.+.309233*t913-.337260*t923+2.51922*t9 
-2.79097*t943+2.16082*t953-.976181*t9*t9 
+.210883*t943*t9-.0169027*t953*t9+7.845538e-6*t9**3) 

5. 3He(n,j3)<: 

Na{(tv) = 6.84713c+8*(l.-.0171094*t912-2.66179*t9+8.27463*t932 
-14.3898*t9*t9+15.6385*t932*t9-10.3337*t9**3 
+3.80177*t912*t9**3-.599790*t9**4-.0139213*t912*t9**4 
+.0140311*t9**5-.00106709*t912*t9**5+1.06709e-6*t9**6) 

6. t{d, n)*He 

Na{(Jv) = 1.78988el2*t9m23*ex(-4.5245/t913) 
/(I. + ((0.129964*t923-0.0482)/(0.5*0.0806))**2) 
*(l.-14.3137899*t913+92.4325675*t923 
-314.645738*t9+641.100355*t943-844.106855*t953 

+752.418564*t9*t9-465.820564*t943*t9 
+202.276143*t953*t9-61.3172473*t9**3 
+12.6913874*t913*t9**3-1.707344*t923*t9**3 
+.134399048*t9**4-.00469341945*t913*t9**4) 

7. 3He(d,p)4He 

Na{(Jv) = 5.67897el2*t9m23*ex(-7.1840/t913) 
/(I. + ((0.206357*t923-0.183)/(0.5*0.256))**2) 
*(l.-8.59410908*t913+31.1979775*t923 

-61.2218616*t9+72.0331037*t943-52.8696341*t953 
+23.7371543*t9**2-5.4569107*t943*t9-.226478266*t953*t9 
+.583380161*t9**3-.190978484*t913*t9**3 
+.031949394*t923*t9**3-.00284146599*t9**4 
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+.106749198e-3*t913*t9**4) 

8. 3He(a,7)^Be 

Na{(Jv) = 3.94207e+6*t9m23*ex(-12.8274/t913) 
*(1. + .185267*t913 - .837432*t923 

+7.23019*t9 - 26.1976*t943 +41.6914*t953 

+19.4465*t9*t9-215.248*t9*t943+422.548*t9*t953 

-412.866*t9**3+176.691*t913*t9**3+45.8891*t923*t9**3 

-100.644*t9**4+54.8984*t913*t9**4-14.1903*t923*t9**4 

+1.48464*t9**5) 

9. t{a,jy Li 

Na{(Jv) = 4.65494351e+6*t9m23*ex(-8.0808/t913) 
*(1. - 12.3956341*t913 + 76.2717899*t923 
-250.678479*t9+446.413119*t943 -289.008201*t953 

-474.786707*t9*t9 +1346.42142*t9*t943 -1503.09444*t9*t953 
+923.138882*t9**3-306.14089*t913*t9**3+42.9886919*t923*t9** 

10. ^Be(n,p)''Li 

Na{(tv) = 5.17900e9*(l.-1.44587*t912+1.12925*t9-.493526*t932 

+.126269*t9*t9-.0194265*t932*t9+.00177188*t9**3 

-.883411e-4*t912*t9**3+.185551e-5*t9**4) 

+4.2994e9*t9m32*ex(-3.713442/t9) 

+1.36949ell*t9m32*ex(-31.332167/t9) 

11. ^Li(j3, a)'*He 

Na{(tv) = 9.19322c8*t9m23*cx(-8.4730/t913) 

*(1. - 2.26222*t913 + 11.3224*t923 

- 27.3071*t9 + 41.1901*t943 - 37.4242*t953 
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+ 18.3941*t9*t9 -3.72281*t9*t943 +2.58125e-2*t953*t9) 
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